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Cosmological perturbations in a healthy extension of Hof ava gravity 
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In Hofava's theory of gravity, Lorentz symmetry is broken in exchange for renormalizability, but 
the original theory has been argued to be plagued with problems associated with a new scalar mode 
stemming from the very breaking of Lorentz symmetry. Recently, Bias, Pujolas, and Sibiryakov have 
proposed a healthy extension of Hof ava gravity, in which the behavior of the scalar mode is improved. 
In this paper, we study scalar modes of cosmological perturbations in extended Hofava gravity. 
The evolution of metric and density perturbations is addressed analytically and numerically. It is 
shown that for vanishing non-adiabatic pressure of matter the large scale evolution of cosmological 
perturbations converges to that described by a single constant, £, which is an analog of a curvature 
perturbation on the uniform-density slicing commonly used in usual gravitational theories. The 
subsequent evolution is thus determined completely by the value of £ . 

PACS numbers: 04.60.-m, 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Recently, Hofava has proposed a new theory of quan- 
tum gravity [l| . The basic idea of the theory is abandon- 
ing the Lorentz invariance, with which the theory is made 
power-counting renormalizable. The broken Lorentz in- 
variance results in a preferred foliation of spacetime by 
three-dimensional spacclike hypersurfaces, and thus al- 
lows to add higher spatial curvature terms to the gravi- 
tational Lagrangian as well as to modify the kinetic term 
of the graviton. Obviously, some notable features of 
Hofava gravity become manifest at high energies, which 
motivates studying the early universe based on Hofava's 
theory @-[3. The effect of the broken Lorentz invari- 
ance, however, persists down to low energies. This not 
only brings interesting consequences regarding in partic- 
ular dark matter in the universe but also causes po- 
tential problems stemming from an additional scalar de- 
gree of freedom of the graviton that inevitably appears 
as a result of the reduced symmetry of the theory 0- 
[Toj . Some of the troubles are cured in the projectable 
version of Hofava gravity 0, and some are still con- 
troversial. Several problems have also been reported at 
a phcnomenological level, including a large isocurvature 
perturbation [l2| and the absence of static stars [l| (see 
also H3). 

It is argued by Bias, Pujolas, and Sibiryakov that origi- 
nal Hofava gravity can be extended in such a way that the 
lapse function N may depend on the spatial coordinate 



(i.e., the theory is not projectable) and terms constructed 
from a 3- vector di In TV are included in the action [l5| , 
which makes N dynamical [ltjj . With the appropriate 
choice of the new terms, the resultant theory will be free 
from the problems and pathologies reported in the liter- 
ature fa-[lfj] . The extended theory could still suffer from 
strong coupling at low energies [17| , but it can be evaded 
by taking into account higher spatial derivative terms in 
the action See also the recent paper for the 

strong-coupling issue in the extended version of Hofava 
gravity. 

So far cosmological perturbations have been explored 
in projectable and non-pro jectable versions of Hofava 
gravity in the vast literature [IH, l2(il - [24| . The aim of this 
paper is to go on in this direction and to study cosmolog- 
ical perturbations in healthfully extended Hofava gravity 
of [15| . We begin with presenting the most general equa- 
tions for scalar perturbations on a flat cosmological back- 
ground in extended Hofava gravity, and then solve them 
both analytically and numerically. As for the other as- 
pects of extended Hofava gravity, spherically symmetric 
solutions have been worked out in [25f . 

The paper is organized as follows. In the next sec- 
tion we review how non-pro jectable Hofava gravity can 
be healthfully extended, following Then, in Sec. Ill, 
the background evolution and perturbation equations are 
given within the framework of extended Hofava gravity. 
In Sec. IV the perturbation equations are solved ana- 
lytically and numerically. We draw our conclusions in 
Sec. V.' 
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II. A HEALTHY EXTENSION OF HORAVA 
GRAVITY 

The constituent variables in Hofava gravity are N , iV,, 
and gij , which correspond to the lapse function, the shift 
vector, and the spatial metric, respectively, in the (3 + 1) 
decomposition of a spacetime. The non-pro jectable ver- 



2 



sion of Hof ava gravity is allowed to include terms con- 
structed from the 3-vector, 



di := dilnN(t,x), 



(1) 



in the Lagrangian. Along with the line of Refs. [15l. |26|. 
we consider the extended Hofava gravity described by 



s = Mk 



J dtd 3 x^/gN (C K -Vfotf.Oi]) 



dtd 3 Xy/g~NC r 



As a result of broken general covariance, a new scalar 
degree of freedom appears in addition to the usual 
helicity-2 polarizations of the graviton, and this scalar 
graviton could be a ghost depending on the value of A. 
Since the structure of the kinetic term ([2]) is the same as 
in the original version of Hofava gravity, the time kinetic 
term of the scalar graviton remains the same. Therefore, 
the condition for avoiding ghosts is 



3A- 1 
A - 1 



> 0. 



(7) 



The "kinetic term" is given by 



(2) 



with the extrinsic curvature defined by Ki 
(<7y - V t N 3 
by V[gij,a,i) 

Vi 
V 2 



VjNi) /2N , while the "potential" is given 
Yjx=i where 
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+r) 2 a i Aa i + r) 3 KVia i + ...), 
V 3 = M Fl \g A RAR + g 5 V t R jk V l R lk 
+ri 4 a i A' 2 a i + n 5 ARVia' + ...) 



(3) 
(4) 
(5) 



Here, the Ricci tensor Rij and Ricci scalar R are con- 
structed from gij, Vj is the covariant derivative with 
respect to gij, A := ViV 4 , and ' := cV In V2 and V3 
we only write explicitly the terms that will be relevant 
to scalar-type, linear perturbations on a flat cosmolog- 
ical background. For instance, terms such as aiajR lJ 
and (ajO 1 ) 2 are not relevant to linear perturbations on 
a flat cosmological background. A parity-violating term 
t % i k RuV jR l k , which appears in the original version of 
Hofava gravity with the detailed valance condition [l[ , is 
not relevant to scalar-type perturbations. The potential 
presented above is therefore the most general one in the 
situation we are considering. C m is the Lagrangian of 
matter. In this paper we simply assume that matter is 
minimally coupled to gravity. 

The symmetry of the theory is invariance under 
the foliation-preserving diffeomorphism transformations: 
t — > t(t), x % — > x l (t,x). Under the infinitesimal trans- 
formation, t — > t + f(t), x % — > x % + £*(t, x), the variables 
transform as 



However, the sound speed squared is negative for srf > 
in original Hofava gravity, indicating that the scalar 
graviton is unstable [§, [22| . This fact itself does not nec- 
essarily mean that the theory is "unhealthy," because 
whether or not an instability really causes a trouble de- 
pends upon its time scale In the above "healthy" 
extension of Hofava gravity, such an instability can be 
made absent from the beginning for stf > by the ap- 
propriate choice of the potential terms of at [la] . 

The Hamiltonian constraint is derived by varying the 
action with respect to N: 



C 



K 



V + 5V 



M 2 1 



0, 



(8) 



where 



SV := 



2 V V i a i -^AV i a i + {£ r AR 



(9) 



We have written here explicitly the terms that are rele- 
vant at zeroth and linear order in our cosmological set- 
ting. The matter energy density is defined as p := 
— £ m — N6C m /SN. Variation with respect to the shift 
vector gives the momentum constraint equations: 



J' 



(10) 



N 
N, 

9ij 



N-fN-fN- fcW, 

Ni - Vi&N, - e j 'v 3 -jvi - e ] gij 

fl, , - 9ijf - .'/-/ V - H,,X,^' 



fNi - fNi 



(0) 



where 7T ij := K ij - XKg^ and J 1 := -N8£ m /8Ni. Fi- 
nally, the evolution equations are derived from variation 
with respect to g^: 



2 (K ik K) - XKKij) - - 9ij C K 



1 d 

W^jdt 
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N 



3 



~AN 9ij - j^ViVjN + /,'., - - 9ij R 



\ (V<V 3 --fl«A)iJ- 



% (V,V, - <^A) Ai? 



pi 



_33_ 

iU pl 



V; V,i2 — A i? 



M- 



pi 



;V,v, //,,Ai 



ff5 

773 \7 « fc 4- 



ViVjAR - A 



^5 

AT 4 



AV fc a fc 



T 

ju pi 



where Tij := £ m gij — 28C m /Sg iJ . Here also we have 
written explicitly only the terms linear in Rij and a,; 
the other possible terms are not relevant in the present 
paper. 

The matter action is invariant under the infinitesimal 
transformation which results in 



d 3 x 



0, 



fl.jT" + Nd t (y/gp) + Nid t (JgJ 1 ) 



(11) 



and 



CLip- 



(Vj-JVi - ViNi) = 0. 



(12) 



Note that the energy conservation equation (fTTj) is given 
by the integration over the whole space, as the gauge 
parameter f{t) is space- independent. 



III. COSMOLOGY OF EXTENDED HORAVA 
GRAVITY 

A. Background equations 

To describe the background evolution of the universe 
we write N = 1, N.i = 0, and = a 2 (t)Sij, assuming 
that the three-dimensional spatial section is flat. The 
Hamiltonian constraint reduces to 



3(3A - 1) 



H 



P 



(13) 



where H := a/ a and p is the background value of the 
energy density of matter, while the evolution equation 
reads 



3A - 1 



3H 2 + 2H 



P 

A/2 ' 



(14) 



where p is the background value of the isotropic pres- 
sure of matter, — <p8? . From the above equa- 
tions we may naturally identify the gravitational constant 
in Friedmann- Robertson- Walker cosmology as 8irG c := 
2Mpf(3A- l)" 1 . 



It follows from Eqs. flT3J and (|T4| that 
p + 3H(p + p) = 0. 



(15) 



Obviously, the energy conservation equation derived 
from (jTTjl . 



J d 3 x a 3 [p + 3H(p + p)] = 0, 



(16) 



is satisfied. Note that in contrast to the case of pro- 
jectable Hofava gravity, we now have the local Hamilto- 
nian constraint and hence "dark matter as an integration 
constant" @ does not appear. 

We see that nothing special happens at the background 
level, except that the gravitational constant G c differs 
from the locally measured value of Newton's constant, 
Gat = il/ P1 2 (l-7 ? /2)- 1 [HI. 



B. Cosmological perturbations 

General scalar perturbations of the lapse function, shift 
vector, and the spatial metric can be written as 

N = 1 + (j>(t,x), Ni = a 2 diP{t,x), 
gij = a 2 [1 - 2^(t, x)} Sij + 2a 2 d l d j E(t, x). (17) 

Note that in the cosmological background of Hofava grav- 
ity we are considering, linear perturbations can be de- 
composed into scalar, vector, and tensor perturbations 
because all the background quantities are homogeneous 
and isotropic. Under a scalar gauge transformation, 
t —> t + f(t), x l — > x l + d l £(t,x), these perturbations 
transform as 



(18) 



By using the spatial gauge transformation we may set 
E = 0. Note, however, that the temporal gauge degree 
of freedom does not help to remove <f> or ip because / is 
a function of t only. 

The linearized Hamiltonian constraint is given by 



3(3A - 1)H hp + H<f>+ iv 2 /3 

V 2 V 2 V 2 
-2—^ + rj—ij,+ —m-\ 



_6p_ 

M 2 X 



0, (19) 



where 



6H(t,x) 



J^V 2 

M 2 , a 2 9 

a 4 



2t?3 V 2 . 
M 2 , a 2 V 

, 2775 V 4 



(20) 
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The momentum constraints at linear order read 

<S.J 



(3\-l)[i> + H<j)) +(A-1)V^ 



M 2 ' 



(21) 



where we write Ji = diSJ. The evolution equations can 
be written as 



{did* 1 A 



(22) 



where 

g T = (3A-1) 



+(\-i)'V[p + 3Hp) + -g TL 



(23) 



and 



Gtl = V 2 ($ + 3Hp) + vv 



- V) + , (24) 



with 



S£(t,x) 



8ff 2 + 3 53 V 2 8^4-355 V 4 
"2"^ ^73 4 V 



Mpj a 2 



A/p, a 4 



(25) 



Neglecting for simplicity the anisotropic stress perturba- 
tion, we have 



Sp 
JU P1 



(26) 
(27) 



Gtl = 0. 

The perturbed energy conservation equation reduces 

to 

J d 3 x a 3 Sp + 3H (Sp + Sp) -3ip(p + p) 

+ J d 3 x a 3 [p + 3H(p + p)] (0 - 3^) - 0. (28) 

The second line vanishes thanks to the background equa- 
tion. Combining the Hamiltonian constraint (|19[) and the 
evolution equations ([25)1 and (|27p. and using the back- 
ground equations, one finds 



Mp" 2 



- z! 

" a 2 
+V 2 



Sp + 3H(8p + 8p)-3ip(p + p) 
2U + HxJj) -r)U + He))) - 



(l-3X)H/3-2H [/3 + 3H/3 



(29) 



Since the right hand side is a total derivative, the first 
line of Eq. ([28)) also vanishes, and hence Eq. ([28)) is auto- 
matically satisfied. Note, however, that the local equa- 
tion ([29)) is stronger than the integrated equation ([28| . 
The momentum conservation equation gives 



5 J + 3HSJ -(p + p)cj> -6p = 0. 



(30) 



This equation can also be derived using the momentum 
constraint ([2"T)) and the evolution equations (J2S)) and ([2"7| 
together with the background equations. Therefore, the 
energy and momentum conservation equations do not 
give rise to any independent equations. 

In order to obtain a closed set of equations, we need 
a matter equation of motion. Suppose that the matter 
equation of motion takes of the form 



Sp + 3H(Sp + Sp) - 3ip(p + p) 



[5J + a 2 (p + p)(3] 



V 2 

— ¥ 5F(t,x), 
or 



(31) 



where the concrete form of ST depends on the matter 
field one is considering. A matter field that respects four- 
dimensional general covariance is conserved locally, and 
consequently, ST = 0. However, in Hofava gravity it is 
natural to consider that matter fields respect only the 
foliation-preserving diffeomorphism, and hence are not 
conserved locally in general. For example, a scalar field 
Lagrangian presented in US HI leads to ST ^ 0. In this 
paper, we assume that SJ 7 arises only from the ultraviolet 
(UV) effect and is suppressed by A/Mj 2 , where M m is a 
typical mass scale. M m is not necessarily of the same 
order of Mp\ and may be different for different matter 
fields, but here we simply assume that the scale M m is 
sufficiently high. Using the momentum constraint (|2ip . 
one obtains 



ST-L + H5H - 2H5£ 



ST_ 

JUpj 



0. 



(32) 



C. Minkowski limit 

Let us check that the perturbation equations derived in 
the previous subsection can reproduce the known result 
in the Minkowski background by taking the limit H — > 0. 
Since we are particularly interested in the stability of a 
scalar graviton in the infrared (IR) regime, we neglect 
the UV terms ST-L and 5£ . Then, the Hamiltonian con- 
straint simply gives <f> = (2/77)^. The trace part of the 
evolution equation ([26)) reduces to s^ip + V 2 $ = 0, while 
the traceless part Q27J| implies $ = (1 — 2/77)1/' • We thus 
obtain 



a/if) - 



n 



^v 2 ^ = 0. 



(33) 



This reproduces the result obtained in [Iff. It can be 
seen that the scalar graviton is stable in the IR regime 
provided that 

< 77 < 2. (34) 
The propagation speed of the scalar graviton is given by 
a _ 2 - 77 1 



77 srf 



(35) 
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In addition to the theoretical bound (JT)) and (|34|) . low 
energy phenomenology can put constraints on the value 
of A and r\. A cosmological bound comes from the fact 
that G c does not coincide with Gn- The bound on the 
difference is derived from the prim ordial abundance of 
He 4 : \G C /G N - 1| < 0.13 [H lH, which in turn gives 
a mild constraint on the parameters. This constraint 
can be further weakened by considering additional rela- 
tivistic degree of freedom and/or large lepton asymmetry 
[29l | to compensate a change in the expansion rate of the 
universe. More stringent constraints will come from the 
parameterized post-Newtonian study. See Ref. [l8j] for a 
preliminary discussion on this point. 



IV. EVOLUTION OF COSMOLOGICAL 
PERTURBATIONS IN THE IR 

In what follows, we are interested in perturbation 
modes whose spatial momenta are much smaller than 
Mpi and M m , so that we ignore any UV terms in gravity 
and matter sectors. 



ip + C* decays also in this case because (1 + e)/2 > p. 
We therefore conclude that the large scale evolution of 
the non-decaying perturbation can be characterized by a 
single constant £*, and 



(40) 



provided that the mode stays in the supcrhorizon regime 
for a sufficiently long time. This indicates that the sub- 
sequent evolution of cosmological perturbations is essen- 
tially determined from the value of . This fact will be 
confirmed by a numerical calculation. 

A remark is in order. In usual gravitational theories 
having general covariance, C is conveniently computed 
from the matter sector by taking the gauge xp = 0. In 
Hofava gravity, however, ip cannot be gauged away be- 
cause the symmetry is reduced and the relevant gauge 
parameter is independent of x. In this sense, two scalar- 
type contributions, a scalar graviton xp and an adiabatic 
matter fluctuation dp, translate to Therefore, we need 
to work out the two physically different contributions in 
order to compute C in Hofava gravity. 



Subhorizon evolution 



A. Superhorizon evolution 

We begin with examining the large scale evolution of 
perturbations analytically by neglecting gradient terms. 
It can be shown from the matter equation of motion with 
5J 7 = that 



C := -V - H 



6p 



(36) 



is conserved on large scales, £ ~ £* = constant, in the 
case of a negligible non-adiabatic pressure perturbation. 
This is the well-known fact when matter is conserved 
locally [27] . Using the Hamiltonian constraint, we obtain 



H 
H 



H<P) ~C*- 



Equation ([32]) on large scales reads 

+ Hep) + cp + H(f> ~ 0, 



(37) 



(38) 



where 28 := 3(A — l)/?/(> 0). The above two equations 
are combined to give 



d 2 y> 

d In a? 



de 
dlr^ 



(39) 



where ip := (tp + C*) exp [i J(l + e)dlna] and e := 
—din H /din a. We may assume e > and de/dlna ~ 
0. (These assumptions are indeed true in a matter- 
dominated universe and a radiation-dominated universe.) 
Then, if (1 — e) 2 — 4e& < 0, %p oscillates and therefore 
i> + C* decays as - a~( 1+£ ) /2. If (1 - e) 2 - 4e,% > 0, $ 
may grow like $ ~ a p with p = — e) 2 — AeSS/2, but 



Let us define 



5 :-- 



Sp + 3H6J 



(41) 



Using the equation of motion of matter, the momentum 
conservation equation, the evolution equations, and the 
Hamiltonian constraint, we find, for a universe domi- 
nated by a fluid with w = p/ p = const., 



5 + (2 - 3w)HS 



w—^-o 

a 1 



+(l + w) 



-w)(l 
3(A-1)H 



+ 3w)H 2 5 
? /V 2 



q>, (42) 



3A- 1 

which coincides with the corresponding equation in gen- 
eral relativity except the last term. The Hamiltonian and 
momentum constraints give 



V 2 



n—4>-2HV 2 p 



_pl_ 

M 2 ' 



(43) 



which can be regarded as a generalized Poisson equation. 
So far no approximation has been made other than ne- 
glecting the UV terms. 

We now focus on the subhorizon evolution of S in a 
matter-dominated universe (w = 0). It is natural to as- 
sume that in the matter-dominated universe the metric 
perturbations are slowly- varying functions of time in the 
sense that their time derivatives give rise to the Hub- 
ble scale: ip t <j> ~ Hip, He/). 1 Under this "quasi-static" 



1 Here we dropped by hand the contribution of a scalar gravita- 
tional wave, but the approximation will be justified by a numer- 
ical calculation. 
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approximation, the evolution equations imply that 

V « 0, P « 0, 



(44) 



on subhorizon scales. Equations (|42[) - (|4"4")) are then ar- 
ranged to give the following equations: 



. V 2 
S + 2HS = -=-<; 



-0 = 4itGnP$, 



(45) 
(46) 



where 8ttGn = (1 — ry/2) -1 . This gravitational con- 
stant coincides with effective Newton's constant defined 
in [l5| (so that it differs from the gravitational constant 
in the Fricdmann equation, G c ). This is a natural exten- 
sion of the result of [15}, where the gravitational field of 
a static point source has been derived. Again, we have 
the relation ij) = <f>, in contrast to the case of Lorentz- 
invariant scalar-tensor theories of gravity. 

Equations (|45|) and (|46|) admit the analytic solution 



(47) 



where £ := G N /G C = (3A - l)/(2 - rj). The first term 
grows in time, and, for £ ^= 1, the growth rate of the 
matter density perturbation is slightly different from the 
standard one (5 ~ t 2 / 3 ). 

Next, let us consider a radiation-dominated universe 
(w = 1/3) and study the evolution of <5 inside the sound 
horizon. From Eq. (|43| we may estimate (V 2 /a 2 )0 ~ 
0(H 2 S), so that we have <5 + H 6 ~ (V 2 /3a 2 )<5 inside the 
sound horizon. This coincides with the corresponding 
equation in general relativity, and the last term hinders 
the growth of the density perturbation in the radiation- 
dominated stage. Moving to the Fourier space and using 
the conformal time defined by dr = dt/a, this equation 
can easily be solved to give S = Cq cos(kT/V3 + 9o), with 
Co and 9q being integration constants. 

Although the behavior of the density perturbation on 
subhorizon scales in the radiation-dominated universe is 
simple as seen above, the metric perturbations evolve 
in a non-trivial manner. Using the generalized Poisson 
equation (|43p. the traceless part of the evolution equa- 
tions (|27j) , and Eq. (|32[) , with some manipulation and the 
approximation k 2 ^> 0(1/t 2 ), we arrive at 



st U" + -(/>' + c 2 a k 2 



6 a p 



7/Af 2 j \k 2 

where ' := d T . The detailed derivation is presented in 
Appendix [A] A solution to this equation is given by 



6fc" 



9 

a z p 



j].st - 3(2 - rj) M, 



[S + 0(S'/k 2 T)] ,(49) 



pi 



where \&(r; k) is a solution to the source- free wave equa- 
tion + (2/r)*' + c 2 fc 2 * = 0, and hence 



cos (cg/cr + 9y) 



(50) 




FIG. 1: Evolution of ip and ^ in a matter-dominated universe. 
Red lines with circles correspond to the initial condition <j>Q = 
ipo{= l), while blue lines with crosses cj>a = — 2ipo(= —2). 
Other initial data are the same and are given by /3o = l X k^ 1 
and xo = 0. The wavenumber is given by k = 0.1xa(to)H(to). 
The parameters are A = 1.05 and rj — 0.1, and a c represents 
the "horizon-crossing time" defined by k — a c H(a c ). 



where Cq, and Qq, are integration constants. The mode 
can be naturally identified as the scalar graviton, and 



it decays as 



inside the horizon, analogously to the 



usual hclicity-2 gravitational wave mode in a cosmologi- 
cal setting. Note that the propagation speed of a scalar 
gravitational wave in a cosmological background is iden- 
tical to c g in the Minkowski background, and its stability 
is ensured for c 2 > 0. The 5 induced part of <f> decays as 
~ a~ 2 . ip can be obtained from the generalized Poisson 



equation (|43[) as 



2 



1 i)s4 - 3(2 + 77) a 2 p 



2k 2 Tj£t - 3(2 



7/0 M 2 ! 



(51) 



The subhorizon evolution in a radiation-dominated uni- 
verse is apparently characterized by four integration con- 
stants, but they are completely determined by a single 
constant, £*, and thus are in fact related to each other. 



C. Numerical solutions 

We have solved numerically the perturbation equations 
in the IR without any other approximations. The pro- 
cedure for doing so is described in some detail in Ap- 
pendix [Bj The results are summarized in Figs. HHU all 
of which reproduce the analytic results obtained above. 

In Fig. Q] we show the evolution of ib and <j> in a matter- 
dominated universe. The initial condition is set by spec- 
ifying four numbers V>o = ip(to), Xo = ip(to), </>o = <£(*o)> 
and Po = (3(to) at some initial time t = to- It can be 
seen that the perturbation evolution starting with the 
two different initial conditions first converges to the so- 
lution (|4T)|) outside the horizon, and then to (|4"4"|) inside 
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FIG. 2: Subhorizon evolution of the density perturbation (di- 
vided by £*) in a matter-dominated universe, compared to the 
analytic solution (|47|l . 
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FIG. 3: Evolution of ip (solid line) and <j) (dashed line) with 
k = 0.01 x a(to)H (to) in a radiation-dominated universe. The 
initial condition is given by ipo — 1 = <f>o, Pa = 0, and xo = 0. 
The parameters are A = 1.05 and r\ = 0.1, and a c is defined 
similarly to the previous figures: k — a c H(a c ). 



the horizon. The evolution of 5 is shown in Fig. [5] for the 
same set of the initial conditions as in Fig. [TJ One can 
confirm that the growth of 5 is in agreement with the 
analytic solution (|T7| . 

The same large scale behavior, ip — ~C* an d — 0, 
is found also in a radiation-dominated universe, as can 
be seen in Fig. [3l The oscillation on small scales in the 
radiation-dominated universe, shown in Fig. 21 is well 
approximated by the analytic solutions presented in the 
previous subsection. The S induced part and VP typically 
have the same amplitude in the plotted region, and \1/ 
dominates for larger r. We have confirmed that different 
initial conditions with the same value of (and the same 
set of the other parameters) result in the identical profile. 
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FIG. 4: Subhorizon evolution of if) and </> in a radiation domi- 
nated universe, compared to the analytic solutions (|49p - (|51|l . 
The initial condition and parameters are the same as in Fig. [3] 

V. CONCLUSIONS 

In this paper, we have studied cosmological perturba- 
tions in a healthfully extended version of Hof ava gravity. 
The most general perturbation equations have been de- 
rived without specifying the matter content. We then 
solved the resultant perturbation equations in the IR 
regime analytically and numerically, for a universe domi- 
nated by a single fluid with vanishing non-adiabatic fluc- 
tuations. It was found that the large scale evolution con- 
verges to that described by a constant £, which is an ana- 
log of the curvature perturbation on uniform density hy- 
persurfaces commonly used in the context of usual gravi- 
tational theories having general covariance. This implies 
that, although the system has two scalar degrees of free- 
dom corresponding to a scalar graviton and an adiabatic 
matter fluctuation, it is sufficient to specify the value of 
C for predicting the late-time evolution of perturbations. 
Our analytic results were confirmed by numerical calcu- 
lations. 

It would be important to revisit the analysis of cos- 
mological perturbations in extended Hofava gravity by 
employing the Hamiltonian formulation in order to get a 
more transparent understanding of the properties of the 
scalar graviton [H, [24[. This issue is left for a further 
study. 

It was pointed out very recently that the IR limit of 
extended Hofava gravity is identical to Einstein-aether 
theory [3(| if the aether vector is restricted to be hy- 
persurface orthogonal [IH, [3l[ . Hypersurface orthogonal 
solutions of Einstein-aether theory are also solutions to 
the IR limit of extended Hofava gravity, though the con- 
verse is not true. This observation is insightful for un- 
derstanding the background cosmological dynamics, as a 
cosmological background aether field is hypersurface or- 
thogonal. For example, the discrepancy between Gn and 
G c is a generic feature of Einstein-aether theory. Note, 
however, that the aether field is no longer hypersurface 
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orthogonal at perturbative order. It would be interesting 
to compare the evolution of cosmological perturbations 
in extended Hofava gravity with that in Einstein-aether 
theory [HIH. 



or equivalently, 



2W 



zlk 2 c 



3 a 2 p 



M. 



(A8) 
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5" + (1 - 3w)-S' + wk 2 S = 1(1 - w)(l + 3w) ( -) S 
a 2 \ a J 



< 1+ ^%~x-i V k ^ (A9) 



Appendix A: Master equations on small scales 

We present a detailed derivation of the master equa- 
tions governing the subhorizon evolution of 5 and <f> for 
general w (=const.). 

Defining convenient variables as u := a 2 ip, v := a 2 0, 
and /3 = a 3 /3, Eq. (|4"3"| can be written as 



a' ~ a p 8 



0. 



(Al) 



while the traceless part of the evolution equations reads 

/3' + v - u = 0, (A2) 

where the conformal time r is used: ' := d T . Equa- 
tion (f3"2")l implies 



a 1 a' 
2—u H v 

a a 



-v = 0. (A3) 



Differentiating Eq. (|A3[) and using Eq. (|A2[) to remove 
/?', we obtain 



2—u' 
a 



+ — (v — u) 



0, 

(A4) 



where we used a" /a, (a 1 /a) 2 k 2 . The expression in the 
square brackets can be written, using Eq. (jAlj) and then 
Eq. (|A2|) . as 



2-v> 
a 



T] k 2 



4 

a p 
2M 2 l 



M 



where 



ka 5 " r a ' S ' 6 



(A5) 
(A6) 



and we used k 2 ^ 0(t~ 2 ) again. Using Eq. 
the expression in the second line can be evaluated as 
~ — £$~ x (a! I a)v' . Thus, we arrive at 



o ' i 2,2 

2 — v + c n k v 
a y 



3 a p 



2 M, 
pi 



(A7) 



We should note that here we do not drop the O ((a'/a) 2 S) 
term because 0(k 2 8) term vanishes for w = 0, in contrast 
to the approximation made in deriving Eqs. (|A4[) . (|A5[) . 
and (|A8[) , where we implicitly assumed that the coeffi- 
cients of Oik 2 ) terms are not too far from the order of 
unity value. In particular, we are not considering the 
case with « 1. Using Eq. (|A9[) . we have 



M = -(1 



a'S'S 
3W) ( «F + 3 



+ 0{</>). (A10) 



Note that S ~ (k/aH) 2 <j> <fi on small scales. 

In a matter-dominated universe, w = 0, the following 
solves Eqs. flXHJ) and (fX9|) : 



6 = C 1 a^ 1+VI + 2 ^^ + C 2 a^ 1 - VT +^/ 4 ,(All) 



k 2 



AnG N p5 + 



(A12) 



where f and Gn are defined in the main text, and ^ is a 
solution to the wave equation '+2(a' / a)*' '+c 2 k 2 ^ = 0, 
and hence can be identified as a scalar gravitational wave. 



decays as 



on small scales. We assumed that 



the amplitude of ^ is small enough not to source S via 
Eq. (fX9|) . Neglecting Eqs. (|Allj) and (|A12j) reproduce 
the result obtained in the main text. 



Appendix B: Solving the perturbation equations 
numerically 



In this appendix we describe the procedure for solving 
the perturbation equations numerically. The UV terms 
SH, 56, and ST a,re all neglected. We use the following 
variables: x '■= 2? : = 8p/(p + p), and J := 5J/(p + p). 
Initial data at t = to arc specified as follows: 

iKto) = ih, x(to) = Xo, 4(to) = <fo, 0(*o) = A>.(Bl) 

Using the constraint equations, we then obtain the initial 
conditions for T> and J at t = to. Choosing to in the 
supcrhorizon regime, we can evaluate from Eq. (|37p . 
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The constraint equations can be solved for 
the Fourier space) as follows: 



and 



= -at 



2^ + (X-l)§^(V + ^HJ) 



we obtain 



V = 3 X -^J- k 2 f3, 
j = 3HclJ + (f> + c 2 s V, 



(B4) 
(B5) 



a 2 H 2 



a?H 2 



rjx + 2Hip + rjHJ 

(3A - 1)§ {V + 2.HJ) 
rl 



(B2) 



(B3) 



where Jz? := rj(k/aH) 2 + 2jz/. These two equations can 
be used to remove (j) and (3 in the other equations, so that 
a system of four first-order equations for -0, x, D, and J 
can be derived. Assuming Sp — c 2 Sp = with c 2 := p/p, 



from the matter equations, and 



st 



X + 'iH X + Hcf>+ (3H 2 - 
— ^ - 4>) = -*tHc 2 s V, 



2H 



(B6) 



from the evolution equations, together with ip = x- In 
the main text the above set of the equations is solved for 
c 2 = w — constant. 
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